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ABSTRACT 

This paper presents a technique for processing accelerometer data 
obtained during the A. E. missions in order to estimate the atmos- 
pheric density profile. A minimum variance, adaptive filter is 
utilized. The trajectory of the probe and probe parameters are in 
a consider mode where their estimates are unimproved but their 
associated uncertainties are permitted an impact on filter behavior. 
Simulations indicate that the technique is effective in estimating 
a density profile to within a few percentage points. 


iii 

preceding page blank not filmed 



Preceding page blank 


CONTENTS 

Page 

ABSTRACT . . . iii 

INTRODUCTION 1 

THE ATMOSPHERIC MODEL 2 

PROCESSING ACCELEROMETER DATA 3 

SIMULATIONS 10 

CONCLUSIONS 12 

ACKNOWLEDGEMENTS 14 

REFERENCES 14 

TABLES 

Table Page 

1 Standard Deviations Used in Simulation 11 


ILLUSTRATIONS 


Figure Page 

1 Percentage Error of Density Estimate versus Altitude 

for First Simulation . 12 

2 Percentage Error of Density Estimate versus Altitude 

for Second Simulation 13 


PRECEDING PAGE BLANK NOT FILMED 


v 



AN ADAPTIVE TECHNIQUE FOR ESTIMATING THE 
ATMOSPHERIC DENSITY PROFILE DURING THE A. E. MISSION 


INTRODUCTION 

The AE-C, D, and E missions are designed to have nominal elliptical orbits with 
perigee altitudes as low as 120 km. At such altitudes the effect of aerodynamic 
heating and atmospheric drag are considerable. This together with the fact that 
the density profile in this portion of the atmosphere is not well known implies 
that data generated by instruments on the probe must be utilized to rapidly and 
accurately estimate atmospheric density as a function of altitude. Since telem- 
etered accelerometer data will be available a certain well precedented procedure 
is suggested. The ideal gas law and the hydrostatic equation can be utilized to 
parameterize the density profile. The accelerometer data can then be processed 
in a differential correction program and the parameters of the atmospheric model 
can be estimated. The accuracy of the technique is limited by the fact that at 
any given time the altitude, velocity, mass, cross sectional area, and drag coef- 
ficient of the probe are all imperfectly known. Also the model used to parame- 
terize the atmospheric density profile introduces further error since it is not a 
completely accurate reflection of reality. The best way to test the accuracy of 
the procedure is to construct a differential correction program for processing 
accelerometer data, generate simulated accelerometer data from an A. E. type 
trajectory using an atmospheric model different from what is assumed in the 
program, add the proper random component to the accelerometer and trajectory 
data, and determine how successful the procedure is in estimating the actual 
density profile. 

This paper is a report on the results of such a test. The atmospheric density 
profile was parameterized by assuming that above a certain reference altitude 
H 0 , temperature is linear with altitude. This assumption together with the ideal 
gas law and the hydrostatic equation permits one to obtain the density profile 
above Ho as a function of the density and temperature at H 0 and the slope of the 
temperature profile above Ho* A recursive, minimum variance filter was de- 
signed to process accelerometer data and estimate these three parameters. The 
velocity and altitude of the probe at any given time along with mass, cross sec- 
tional area, and drag coefficient are treated in the filter as so called "consider" 
parameters. Thus the uncertainties of these parameters were permitted to have 
an impact on filter performance but their estimates were left unimproved. An 
adaptive feature was utilized to dynamically weight the accelerometer observa- 
tions as a function of residuals. The addition of this feature was found to sub- 
stantially improve the filter’s performance. 
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A numerical integration scheme was used to generate A. E. type trajectories 
and the correct values of altitude, scalar velocity, and drag deceleration were 
recorded for various time points. With the aid of a random number generator, 
appropriate random components were added to these values as well as to probe 
mass, cross sectional area, and drag coefficient. The resultant simulated data 
was then processed by the filter and the estimated density profile was compared 
to the correct density profile. The process was repeated for two different at- 
mospheres. The first atmosphere possessed a linear temperature profile above 
H 0 as assumed by the filter atmospheric model. The results of this simulation 
demonstrate the optimum accuracy to be expected from the procedure. The 
second simulation was performed with an atmospheric model whose temperature 
profile above H 0 was highly non linear. This simulation provides the results to 
be expected in the likely case that significant modeling errors are present. 

Succeeding sections provide the mathematics of the atmospheric model and the 
recursive filter followed by details of the simulations. 


THE ATMOSPHERIC MODEL 

Assume that above a reference altitude H 0 , the Earth's atmosphere is spherically 
symmetric and perfectly mixed. Then the ideal gas law permits us to write 

P(h) = p(h)RT(h)/m (1) 


where P (h), p (h), and T (h) are respectively hydrostatic pressure, temperature, 
and density at altitude h. The symbols R and m respectively represent the ideal 
gas constant and the average molecular mass. The hydrostatic equation provides 


dP (h) 
dh 


- g P (h) 


( 2 ) 


where g is the gravitational constant (assumed independent of h). From Equa- 
tions 1 and 2, one can quickly derive 


1 dp (h) „ , f gm dT (h) 

p (h) dh L R dh J 


(3) 


Differential Equation 3 can be solved for p (h), h > H 0 , only if p 0 = p (H 0 ) is 
known and if T (h), h > H 0 is given. It is assumed in this paper that T (h) is 
linear above H 0 with T 0 = T (H 0 ) and S the slope of T (h). The parameters p 0 , 
T 0 and S serve to define p (h), h > H 0 . To see how, integrate Equation 3 from 
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H 0 to some higher altitude h. This immediately yields 


L °g e 


Pi h) 

. Pq - 


= - Log e 


T (h) 



( 4 ) 


and by taking advantage of the linear form of T (h) one can easily derive from 
Equation 4, 

l +gm/R s 

, h > H 0 (5) 


Pi h) = P 0 


! T 0+ S(h-H 0 ) 


PROCESSING ACCELEROMETER DATA 

The scalar acceleration due to drag experienced by a probe in an atmosphere is 
given by 


a - C V 2 p (h) (6) 

where p (h) is the atmospheric density at altitude h, V is the scalar velocity of 
the probe and C is given by 


C =-~-A C D /n (7) 

where A is the cross sectional area of the probe, C D is the drag coefficient and 
n is the probe mass. In the case of an A. E. probe, the drag coefficient is de- 
termined analytically as a function of angle of attack by assuming the probe to 
be a cylinder. The cross sectional area and mass are continuously monitored. 
A tri-axial set of accelerometers are assumed to be installed on the probe. 
These accelerometers measure three mutually perpendicular components of the 
acceleration due to drag. The root mean square value of the three accelerom- 
eters provides an observation of the scalar acceleration due to drag as symbol- 
ized on the left side of Equation 6. Combining Equations 5 and 6 yields 

a - CV 2 p Q 

The altitude h and scalar velocity V of the probe at the instant an accelerometer 
reading is obtained is assumed available as the output of an orbit determination 
program. The problem can now be stated. Given a series of measurements of 
the left side of Equation 8, determine the most efficient filter for processing 


1 +gm/R S 


L 0 


T n +S(h-H 0 ) 


( 8 ) 
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these observations to obtain estimates of p 0 , T 0 , and S. The development of 
such a filter must take into account uncertainties in C, h, and V as well as un- 
certainties in the observations. For reasons of computational convenience it is 
desired that the resultant filter be both linear and recursive. In order to utilize 
results of linear filtering theory it is necessary to assume a linear approxima- 
tion to Equation 8. Let Tj be the time at which the i th accelerometer measure- 
ment is taken. Then define the following symbols 


r p o i 


[hi I 

T 

. Y i = 

V i 

S 


L c i J 


(9) 


where hj, V, , and Cj are respectively the values of h, V, and C of Equation 8 
at time Tj . Equation 8 can be symbolized as 


Oj=/(X,Yj) (10) 

where Oj is the true acceleration due to drag at time Tj . Let X n and Yj ;I1 be 
nominal estimates of X and Yj which may be thought of as the best estimates of 
X and Yj immediately before the i th accelerometer measurement is processed. 
A linear approximation to Equation 10 must take the form 


a { = Gj X + G 2 Yj + G 3 


( 11 ) 


For our purposes it is not necessary to know a precise value for G 3 . The esti- 
mates used for matrices Gj and G 2 are 


G i = 


a/(x, Yj) 


9 X 


x = x n ,Y.= Y. n 

n i i,n 


G 2 = 


9 f(X, Yj) 


9 Y: 


( 12 ) 


X = X n > Y r Y i,n 


At time Tj , an estimate Yj of Yj is available with statistics. 


Yj = Yj+Tj.Ety) = 0,E(7 i 7 1 T ) = P i ,E(7 i 7 J T ) = 0, i * j (13) 

A 

Since the filter is assumed recursive, an estimate Xj_j of X based on observa- 
tions obtained previous to time Tj is available with statistics 


Xj.i = X + a i _ 1 , E (a i _ 1 ) = 0, E (a j _ 1 , aj _ l ) = P M 


(14) 


The matrices pj and P i _ 1 are presumed known. At time Tj a direct observation 
Oj of cij is obtained with statistics 


ttj = a i + v i + tj, E (vj) = E (tj) 

E(^ 2 ) = 


= 0, E (i>j Tj) = E (i^j v.) = E (tj Tj) = 0, i^j 
Q> E (t 2 ) = Rj (15) 
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The random variable v x represents the inherent instrument noise and its variance 
Q is considered as known. The random variable r, represents observation noise 
due to modeling errors such as inaccuracies in Equation 8 and errors caused by 
the linearization symbolized by Equation 11. Its variance Rj is not known 
a priori but is determined by the residuals of the filter in a way which will be 
discussed later. 

Since the filter is constrained to be linear, the estimate X, is obtained by pro- 
cessing observation cq must take the form 

\ - H ! a, + H 2 Xj., + H, Y, + H 4 (16) 

where H lt H 2 , H 3 , and H 4 are matrices of the proper dimension. At this point 
the unbiasedness condition is imposed on the filter. In short, we insist that 

E(Xj) = X 07) 

Condition 17 has important consequences. By using Equations 11 through 15, 
Equation 16 can be written as 

(18) 

X, = [Hj G t +H 2 ] X + [Hj G, + H 3 ] Yj + Hj G 3 + H 4 + H, v. + H 2 + H 3 7j 

The only way to insure that Equation 17 is satisfied is to impose conditions 

Hj G t + H 2 = I, H t G 2 + H 3 = 0, Hj G 3 + H 4 = 0 (19) 

Equations 19 permits us to represent matrices H 2 , H 3 , and H 4 in terms of Hj. 
Equation 16 can then be written as 

X( = X,., +H, [Oj-G, Xj., -G 2 Yi-G,! (20) 

A A 

The expression G t Xj_j + G 2 Yj + G 3 may be viewed as the best estimate of cq 
previous to processing oq and assuming linear approximation 11 to Equation 10. 
Thus Equation 20 represents the filter estimate as the sum of the previous 
estimate Xj_ , and the product of the vector Hj by the difference between the ob- 
served value a; of cq and the computed best estimate value of oq based on esti- 
mates Xj_, and Yj of X and Y, . This suggests that we define 

a? = / (X i _ 1 , Y) (21) 

where the function is defined by Equation 10 and is explicitly represented by 
the right side of Equation 8. Equation 20 can be restated as 

Xj = Xj.i +H i (“i- a i C) (22) 
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where H, is an as yet undetermined vector. To determine H, we impose the so 
called minimum variance condition on Xj . Define 


p i = 


^(Xj - X) (Xj - X) 



(23) 


A A 

Since the expected value of Xj is X, Pj represents the covariance matrix of X, . 
We intend to choose a value for H, which minimizes the trace of Pj. The result, 
ignoring the non linearity of Equation 21, will be a linear, unbiased, minimum 
variance filter. Equations 18 and 19 permit one to write 


Xj - X+H t t’j + (1 - Hj G, ) a i _ ] -H, G 2 7j + H, Tj 


(24) 


Hence 


Pj = Hj Q H; + Rj Hf + (I - Hj Gj ) R., (I - H t )' + H, G 2 p. G 2 H} (25) 


Define 


5 . =Q+Rj + GiPj_ 1 G[ + G 2 P i G] 


5 2 = P i-i 


(26) 


Then Equation 25 can be written as 


P i = § I H 1 +P i-1 


(27) 


The problem now is to obtain the value of H j which minimizes the trace of Pj. 
Assume the following change of variables 


Hl = m; 1+a 


(28) 


Then Equation 27 can be written 

P i = P i-1 “ 5 2 § I 1 5 I +A5 1 A T < 29 ) 

Since 5j is necessarily a positive number, the trace of Pj is minimized when 
A - 0 and the value of Hj which provides a minimum variance filter is 


( 30 ) 



And from Equations 29 and 30, the covariance matrix of the minimum variance 
filter is 


P. 

i 


P i-, -«.■ 



(31) 


The linear, unbiased, minimum variance, recursive filter can be given as 


• .<=/(X i . 1 ,Y i ) 

(32a) 

Q + R i + G , p i-, G I + G 2^i G I 

(32b) 

6 2 = P i-. G I 

(32c) 

H = 6 2 5' 1 

(32d) 

Xj = Xj., +H(a i -a^) 

(32e) 

P i = P i-.- H5 I 

(320 


Equations 32 are not complete since the value of Rj is undetermined. A mode of 
estimating Rj as a function residuals will now be developed. The procedure will 
again be a recursive one. Hence we assume that an estimate Rj_j of Rj_, was 
obtained at time T_, with statistics 

Rj., = Rj_, + *j_, , E («,_, ) = 0,E(^) = j3j_, (33) 

Assume further that the relationship between Rj and R is 

Rj = R j _ 1 +n { , E (7Tj) = 0 ,E(tt?) = co (34) 

Define the residual of the filter at time Tj as 


Using the linear approximation implied by Equation 11, Equation 35 can be 
written 


e i = I 'i +r i- G l “i-l - G 2^i 

Next define Xj as the variance of e t 

\ = E ( e i 2 ) 


(36) 


(37) 
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Equations 36 and 37 lead to 

Xj = Q + Rj + G t Pj j Gj + G 2 p. G* (38) 

Equations 34 and 38 can be viewed respectively as linear state and observation 
equations. If an unbiased observation of Xj could be obtained it could then be 
processed in a standard minimum variance fashion and a linear, unbiased, mini- 
mum variance, recursive filter could be obtained for estimating Rj . It can be 
shown that the square of the residual ej obtained at time Tj is a chi-square ran- 
dom variable with expectation equal to \ and variance equal to 2Xj 2 . See [1] for 
details. The number e\ can now be treated as an observation of Xj. Its variance 
can be estimated from Equation 38 by substituting Rj.jfor Rj on the right side. 
The result is 


E [(e? ) (e 2 )] * 2 [Q+R M +Gj G[+G 2 PjG^] 2 (39) 

The resultant estimator for Rj is given by the following recursion relations 

X = Q + Rj.j + Gj P M G] + G 2 p. G] (40a) 

q = 2 x 2 (40b) 

z = (3j_j + co (40c) 

K = z (z + q)" 1 (40d) 

Rj = R M +Jf(e 2 -X) (40e) 

Pi = zd-W) 

Again the details can be found in [ 1] . It is instructive to see how Equations 40 


interface with Equations 32 to form the complete filter. Given the mathematical 


model of Equations 11 through 15 and given Xj_j , Yj_j , Pj_j , 
minimum variance recursive filter is 

R i_i , P j_j , the 

=/ ( X i_i> Y i) 

(41a) 

A C 

6- = a. - a. 

i i i 

(41b) 

X = Q + Rj_j + Gj Pj_j Gj + G 2 p. G^ 

(41c) 

q = 2x 2 
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Z = + CO 

(4 1 e) 


X = z (z + q) -1 

(410 

A 

R. 

i 

= R m + X (e\ - X ) 

(4 lg) 


n 

N 

/— s 

1 

*5 

(4 1 h) 

Q + 

\ + G 1 P i-1 G 1 + G 2 ^ G I 

(4 1 i) 


5 2 = P i-, G [ 

(4 1 j) 


h = s 2 s; 1 

(41k) 


= Xj_j +H(a i -a c ) 

(411) 


P i = P i-,- H5 I 

(41m) 


The structure of this set of recursion relations may become clearer if one views 
Equations 41 as a pair of interlocking minimum variance filters. The inner 
filter defined by Equations c through h estimates observation noise variance R, 
and the outer filter defined by Equations a through b and i through m estimates 
the state X. The outer filter passes the residual to the inner filter which 
processes ef as an observation and estimates Rj . The number R, is then treated 
by the outer filter as observation noise variance and an estimate X } of X is 
obtained. 


The filter is adaptive in the sense that a few large residuals will cause the inner 
filter to return a large estimate of observation noise variance with a resultant 
downgrading of the impact of observations on the estimate. In this sense the 
filter adapts its behavior to the presence of modeling errors. 


An adaptive feature when used in this way amounts to a dynamic weighting scheme 
which automatically downgrades the impact of data points which are associated 
with large residuals. Thus the necessity of editing the data before processing is 
eliminated. This is important with reference to real time applications. Notice 
also that the filter defined by Equations 41 does not utilize the data to produce an 
improvement in the estimate Yj of Y, . This means that the accelerometer data 
is used to improve the estimate of the atmosphere rather than the actual trajec- 
tory or the parameters associated with the probe. The uncertainties of the tra- 
jectory and the probe parameters as given by covariance matrix p { do, however, 
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influence the estimate Xj of X. The parameters which define the vector Yj are 
said to be in a "consider mode" rather than a "solved for" mode. 


SIMULATIONS 

Although the results of linear filtering theory were used to generate the results 
of the previous section, the filter defined by recursion relations 41 is non linear. 
Consequently it is not possible to evaluate its performance by ensemble studies. 

If one is to know how such a filter will perform during the actual mission, simu- 
lations must be performed. For the first simulation to be discussed a numerical 
integrator was employed to generate an A. E. type trajectory. The atmospheric 
model was obtained by assuming the ideal gas law and the hydrostatic equation 
and a temperature profile which is a polygonal arc with four breakpoints the 
highest of which at the reference altitude. Since for this model the temperature 
above the reference altitude is linear with altitude, no modeling errors were 
present for this simulation. At discrete time points, so called true values of 
altitude, velocity, and acceleration due to drag were recorded. A random num- 
ber generator was utilized to add appropriate random components to these num- 
bers. Random components were also added to the values of the cross sectional 
area, mass, and drag coefficient that were used to generate the true data. Re- 
cursion relations 41 were translated into a computer program which processed 
the simulated data to yield an estimate of the atmospheric density profile. The 
performance of the filter was measured by the percentage error between the 
estimated density and the actual density as a function of altitude. 


The perigee altitude for the trajectory was 135 km. From 7 minutes before 
perigee until 7 minutes after perigee, a total of 175 accelerometer measure- 
ments were obtained. The uncertainty figures used in the simulations are pro- 
vided in Table 1. The reference altitude was chosen at 100 km. For this esti- 
mation technique it is necessary to begin with an a priori estimate of the 
atmospheric density profile. For this simulation a very poor a priori estimate 
was used. This was done for two reasons. First, the filter's performance 
would not be very impressive if in order to obtain a good estimate of the density 
profile it must begin with a good estimate. The second reason is that in pseudo- 
real time applications such as this the question of the stability of the filter be- 
comes important. If a poor a priori estimate is provided to a recursive fil ter 
a frequent result is filter divergence and a new a priori estimate must be tried. 
This involves a searching procedure which can become lengthy. 


The results of the simulation are displayed in Figure 1 which is a plot of the 
estimated density minus the true density divided by the true density and multiplied 
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Table 1 


Standard Deviations Used in Simulation 


Parameter 

Standard Deviation 

Altitude 

500 m 

Velocity 

4 m/sec 

Accelerometer Reading 

5 (10) -5 m/sec 2 

Probe Cross Sectional Area 

0.16 m 2 

Probe Mass 

6.5 kg 

Probe Drag Coefficient 

0.13 


by 100. It is seen that between 130 km and 200 km the filter succeeds in esti- 
mating the density with an average error of less than 2%. In the same region the 
a priori estimate of the density profile would have introduced an average error of 
well over 600%. 


In the simulation which produced Figure 1 the assumption that above the refer- 
ence altitude temperature is linear is satisfied. During the actual mission this 
is not likely to be the case. In fact present information [2 ] indicates that tem- 
perature above 100 km is likely to possess a negative second derivative with 
regard to altitude. This difficulty could be solved by modeling several break- 
points above 100 km and solving for the extra slopes. But this leads to a more 
complicated filter. It is interesting to observe how the present filter performs 
in estimating the density profile of an atmosphere whose temperature profile is 
not linear above the reference altitude. The second simulation discussed in this 
paper was performed identically to the first with the exception that the atmos- 
pheric model had six temperature breakpoints the last two of which were at 
125 km and 160 km. The slopes were chosen so that the resulting temperature 
profile resembled the curves shown in [2 ]. The results are shown in Figure 2 
which can be interpreted in the same way as Figure 1. The average error be- 
tween 130 km and 200 km has increased to about 3%. This suggests that an 
atmospheric model more elaborate than what is used in the present filter may 
not be necessary in order to adequately estimate the density profile. The same 
simulations were completed without the adaptive feature. The average error 
between 130 km and 200 km for simulation one was about 3% and for simulation 
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Figure 1. Percentage Error of Density Estimate versus 
Altitude for First Simulation 


2 about 4.5%. The adaptive feature appears to be well worth the added complex- 
ity it introduces. This added complexity is not great as seen by the fact that 
with the adaptive feature included the filter is capable of processing over 1,500 
accelerometer readings per minute on a modern high speed computer. 


CONCLUSIONS 

Accelerometer data collected during the A. E. mission along with trajectory 
data can be processed in order to yield an estimate of the atmospheric density 
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Figure 2. Percentage Error of Density Estimate versus 
Altitude for Second Simulation 


profile that should be accurate to within a few percentage points. In this paper 
the filter used to process accelerometer data was minimum variance, recursive, 
and with an adaptive observation noise estimator. The altitude, velocity, and 
probe parameters were in a consider mode where their estimates were unim- 
proved but their associated uncertainties were permitted to have an impact on the 
filter estimate. The atmospheric model used by the filter characterize the 
density profile with just three parameters. This appears to be adequate since >' 
simulations show that even when the modeling assumptions are seriously violated 
the filter performs well. Simulations also show that the filter converges to a 
good estimate even when the a priori estimate of the density profile is in error 
by several hundred percentage points. The adaptive estimation of observation 
was found to add substantially to the accuracy of the technique. 


As designed the filter should be relatively insensitive to data dropouts and to 
badly biased data points. Its speed should be adequate for real time applications 
and its accuracy may be adequate for post flight analysis. 
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